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1.1 Introduction 

The story of this metric begins with a paper by Alexei Zinovievich Petrov 
(1954) where the simultaneous invariants and canonical forms for the 
metric and conformal tensor are calculated at a general point in an Ein- 
stein space. This paper took a while to be appreciated in the West, 
probably because the Kazan State University journal was not readily 
available, but Felix Pirani (1957) used it as the foundation of an arti- 
cle on gravitational radiation theory. He analyzed gravitational shock 
waves, calculated the possible jumps in the Riemann tensor across the 
wave fronts, and related these to the Petrov types. 

I was a graduate student at Cambridge, from 1955 to 1958. In my last 
year I was invited to attend the relativity seminars at Kings College in 
London, including one by Felix Pirani on his 1957 paper. At the time 
I thought that he was stretching when he proposed that radiation was 
type N, and I said so, a rather stupid thing for a graduate student with 
no real supervisor to dcQ. It seemed obvious that a superposition of type 
N solutions would not itself be type N, and that gravitational waves near 
a macroscopic body would be of general type, not Type N. 

Perhaps I did Felix an injustice. His conclusions may have been over- 
simplified but his paper had some very positive consequences. Andrzej 
Trautman computed the asymptotic properties of the Weyl tensor for 
outgoing radiation by generalizing Sommerfeld's work on electromag- 
netic radiation, confirming that the far field is Type N. Bondi, M.G.J, 
van der Burg and Metzner (1962) then introduced appropriate null coor- 



f My nominal supervisor was a particle physicist who had no interest in general 
relativity. 
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dinates to study gravitational radiation in the far zone and related this 
to the results of Petrov and Pirani. 

In 1958 I went to Syracuse University as a research associate of Peter 
Bergman. While there I was invited to join Joshua Goldberg at the 
Aeronautical Research Laboratory in Dayton Ohic0. There was another 
relativist at the lab, Dr Joseph Schell, who had studied Einstein's unified 
field theory under Vaclav Hlavaty. Josh was about to go on study leave 
to Europe for a few months, and did not want to leave Joe by himself. 

Before he left, Josh and I became interested in the new methods that 
were entering general relativity from differential geometry at that time. 
We did not have a copy of Petrov's paper so our first project was to re- 
derive his classification using projective geometry, something which was 
being done by many other people throughout the world at that time. In 
each empty Einstein space, (£, the conformal tensor determines four null 
"eigenvectors" at each point. The metric is called algebraically special 



if two of these eigenvectors coincide. This vector is then called a 



principal null vector (PNV) and the field of these is called a principal 
null "congruence". 

After this, we used a tetrad formulation to study vacuum Einstein 
spaces with degenerate holonomy groups (Goldberg and Kerr (1961), 
Kerr and Goldberg (1961)). The tetrad used consisted of two null vectors 
and two real orthogonal space-like vectors, 



We proved that the holonomy group must be an even dimensional sub- 
group of the Lorentz group at each point, and that if its dimension is 
less than six then coordinates can be chosen so that the metric has the 
following form: 



f There is a claim spread on internet that we were employed to develop an antigravity 
engine to power spaceships. This is rubbish! The main reason why the The US 
Air Force had created a General Relativity section was probably to show the 
navy that they could also do pure research. The only real use that the USAF 
made of us was when some crackpot sent them a proposal for antigravity or for 
converting rotary motion inside a spaceship to a translational driving system. 
These proposals typically used Newton's equations to prove non-conservation of 
momentum for some classical system. 

X The term algebraically degenerate is sometimes used instead. 




ds 2 = (cj 1 ) 2 + (lu 2 ) 2 + 2^uj\ 



ds 2 = dx 2 + dy 2 + 2du(dv + pdx + ^(uj — p :X v)du), 
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where both p and uj are independent of and 

P,xx ~\~ P .yy 
U,xx + LO.yy = %P,ux ~ %PP,xx ~ (P,x) 2 + (P,y) 2 

This coordinate system was not quite uniquely defined. If p is bilinear in 
x and y then it can be transformed to zero, giving the well-known plane- 
fronted wave solutions. These are type N, and have a two-dimensional 
holonomy groups. The more general metrics are type III with four- 
dimensional holonomy groups. 

In September 1961 Joshua joined Hermann Bondi, Andrzej Traut- 
man, Ray Sachs and others at King's College in London. By this time 
it was well known that all such AS spaces possess a null congruence 
whose vectors are both geodesic and shearfree. These are the degenerate 
"eigenvectors" of the conformal tensor at each point. Andrzej suggested 
to Josh and Ray how they might prove the converse. This led to the 
celebrated Goldberg-Sachs theorem (see Goldberg and Sachs (1962)): 

Theorem 1 A vacuum metric is algebraically special if and only if it 
contains a geodesic and shearfree null congruence. 

Either properties of the congruence, being geodesic and shear-free, or a 
property of the conformal tensor, algebraically degeneracy, could be con- 
sidered fundamental with the others following from the Goldberg-Sachs 
theorem. It is likely that most thought that the algebra was funda- 
mental, but I believe that Ivor Robinson and Andrzej Trautman (1962) 
were correct when they emphasized the properties of the congruence in- 
stead. They showed that for any Einstein space with a shear-free null 
congruence which is also hypersurface orthogonal there are coordinates 
for which 

ds 2 = 2r 2 p- 2 d(d(-2dudr- (AlnP - 2r(lnP), u - 2m(u)/r)du 2 , 
where £ is a complex coordinate, — (x + iy)/\^2, say, so that 

2d(d( = dx 2 + dy 2 . 
The one remaining field equation is, 

AA(lnP) + 12m(lnP) iU - 4m,„ = 0, A = 2P 2 d c d^. 
The PN\{±| is k = k^d^ = d r , where r is an affine parameter along 

§ The simple way that the coordinate v appears was to prove typical of these alge- 
braically special metrics, 
f The letter k will be used throughout this article to denote the PNV. 



4 R. P. Kerr 




Fig. 1.1. Ivor Robinson and Andrzej Trautman constructed all Einstein spaces 
possessing a hypersurface orthogonal shearfree congruence. Whereas Bondi 
and his colleagues were looking at spaces with these properties asymptot- 
ically, far from any sources, Robinson and Trautman went a step further, 
constructing exact solutions. (Images courtesy of Andrzej Trautman and the 
photographer, Marek Holzman) 

the rays. The corresponding differential form is k = k^dx^ = du, so 
that k is the normal to the surfaces of constant u. The coordinate u 
is a retarded time, the surfaces of constant r, u are distorted spheres 
with metric ds 2 — 2r 2 P~ 2 d(d( and the parameter m(u) is loosely con- 
nected with the system's mass. This gives the complete solution to the 
Robinson-Trautman probleirjjl 

In 1962 Goldberg and myself attended a month-long meeting in Santa 
Barbara. It was designed to get mathematicians and relativists talking 
to each other. Perhaps the physicists learned a lot about more mod- 
ern mathematical techniques, but I doubt that the geometers learned 
much from the relativists. All that aside, I met Alfred Schild at this 

\ In the study of exact solutions, "solving" a problem usually means introducing 
a useful coordinate system, solving the easier Einstein equations and replacing 
the ten components of the metric tensor with a smaller number of functions, 
preferably of less than four variables. These will then have to satisfy a residual set 
of differential equations, the harder ones, which usually have no known complete 
solution. For example, the remaining field equation for the Robinson-Trautman 
metrics is highly nonlinear and has no general solution. 
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Fig. 1.2. Ezra T Newman, with T. Unti and L.A. Tambourino, studied the 
field equations for diverging and rotating algebraically special Einstein spaces. 

conference. He had just persuaded the Texas state legislators to finance 
a Center for Relativity at the University of Texas, and had arranged 
for an outstanding group of relativists to join. These included Roger 
Penrose and Ray Sachs, but neither could come immediately and so I 
was invited to visit for the 62-63 academic year. 

After Santa Barbara, Goldberg and myself flew to a conference held 
at Jablonna near Warsaw. This was the third precursor to the trien- 
nial meetings of the GRG society and therefore it could be called GR3. 
Robinson and Trautman (1964) presented a paper on "Exact Degenerate 
Solutions" at this conference. They spoke about their well-known solu- 
tion and also showed that when the rays are not hypersurface orthogonal 
coordinates can be chosen so that 

ds 2 = -P 2 [(d£ - ak) 2 + (dr) - bk) 2 ] + 2dpk + ck 2 , 

where, as usual, k is the PNV. Its components, k a , are independent of 
p, but a, 6, c and P may be functions of all four coordinates. 

I was playing around with the structure of the Einstein equations 
during 1962, using the new (to physicists) methods of tetrads and differ- 
ential forms. I had written out the equations for the curvature using a 
complex null tetrad and self-dual bivectors, and then studied their inte- 
grability conditions. In particular, I was interested in the same problem 
that Robinson and Trautman were investigating where k was not a gra- 
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dient, i.e twisting, but there was a major road block in my way. Alan 
Thompson had also come to Austin that year and was also interested 
in these methods. Although there seemed to be no reason why there 
should not be many algebraically special spaces, Alan kept quoting a 
result from a preprint of a paper by Newman, Tambourino and Unti 
(1963) in which they had "proved" that the only possible space with a 
diverging and rotating PNV is NUT space, a one parameter generaliza- 
tion of the Schwarzschild metric. They derived this result using the new 
Newman-Penrose spinor formalism (N-P). Their equations were essen- 
tially the same as those obtained by people such as myself using self-dual 
bivectors: only the names are different. I could not understand how the 
equations that I was studying could possibly lead to their claimed result, 
but could only presume it must be so since I did not have a copy of their 
paper. 

In the spring of 1963 Alan obtained a preprint of this paper and loaned 
it to me. I thumbed through it quickly, trying to see where their hunt 
for solutions had died. The N-P formalism assigns a different Greek 
letter to each component of the connection, so I did not try to read it 
carefully, just rushed ahead until I found what appeared to be the key 
equation, 



where the rij were all small integers, between -4 and +4. Their sum was 
not zero so this gave a = 0. I had no idea what a represented, but its 
vanishing seemed to be disastrous and so I looked more carefully to see 
where this equation was coming from. Three of the previous equations, 
each involving first derivatives of some of the connection components, 
had been differentiated and then added together. All the second deriva- 
tives cancelled identically and most of the other terms were eliminated 
using other N-P equations, leaving (|1.1[) . 

The mistake that Newman et al. made was that they did not notice 
that they were simply recalculating one component of the Bianchi iden- 
tities by adding together the appropriate derivatives of three of their 
curvature equations, and then simplifying the result by using some of 
their other equations, undifferentiated. The result should have agreed 
with one of their derived Bianchi identities involving derivatives of the 
components of the conformal tensor, the ^ functions, giving 



i(ni + n- 2 + n 3 )a 2 = 0, 



(1.1) 



m + n-2 + n 3 = 0. 



(1.2) 



In effect, they rediscovered one component of the identities, but with 



The Kerr and Kerr-Schild metrics 



7 



numerical errors. The real fault was the way the N-P formalism confuses 
the Bianchi identities with the derived equations involving derivatives 
of the variables. 

Alan Thompson and myself were living in adjoining apartments, so I 
dashed next door and told him that their result was incorrect. Although 
it was unnecessary, we recalculated the first of the three terms, ni, 
obtained a different result to the one in the preprint, and verified that 
(|1.2[) was now satisfied. Once this blockage was out of the way, I was then 
able to continue with what I had been doing and derive the metric and 
field equations for twisting algebraically special spaces. The coordinates 
I constructed turned out to be essentially the same as the ones given in 
Robinson and Trautman (1962). This shows that they are the "natural" 
coordinates for this problem since the methods used by them were very 
different to those used by me. Ivor loathed the use of such things as N-P 
or rotation coefficients, and Andrzej and he had a nice way of proving 
the existence of their canonical complex coordinates £ and C- I found 
this same result from one of the Cartan equations, as will be shown in 
the next section, but I have no doubt that their method is more elegant. 
Ivor explained it to me on more than one occasion, but unfortunately I 
never understood what he was sayin^! 

At this point I presented the results at a monthly Relativity conference 
held at the Steven's Institute in Hoboken, N.J. When I told Ted Newman 
that Ijl.ljl should have been identically zero, he said that they knew that 
ri\ was incorrect, but that the value for ni given in the preprint was a 
misprint and so (|1.2p was still not satisfied. I replied that since the sum 
had to be zero the final term, must also be incorrect. Alan and I 
recalculated it that evening, confirming that (|1.2[) was satisfied 



1.2 Discovery of the Kerr metric 

When I realized that the attempt by Newman et al. to find all rotating 
AS spaces had foundered and that Robinson and Trautman appeared to 
have stopped with the static ones, I rushed headlong into the search for 
these metrics. 

Why was the problem so interesting to me? Schwarzschild, by far 
the most significant physical solution known at that time, has an event 

f While writing this article I read their 1962 paper and finally understood how they 

derived their coordinates. It only took me 45 years. 
| Robinson and Trautman (1962) also doubted the original claim by Newman et al. 

since they knew that the linearized equations had many solutions. 
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horizon. A spherically symmetric star that collapses inside this is forever 
lost to us, but it was not known whether angular momentum could stop 
this collapse to a black hole. Unfortunately, there was no known metric 
for a rotating star. Schwarzschild was an example of the Robinson- 
Trautman metrics, none of which could contain a rotating source as 
they were all hypersurface orthogonal. Many had tried to solve the 
Einstein equations assuming a stationary and axially symmetric metric, 
but none had succeeded in finding any physically significant rotating 
solutions. The equations for such metrics are complicated nonlinear 
PDEs in two variables. What was needed was some extra condition that 
would reduce these to ODEs, and this might be the assumption that the 
metric is AS. 

The notation used in the rest of this paper is fairly standard. There 
were two competing formalisms being used around 1960, complex tetrads 
and spinors. I used the formejjj], Newman et al. the latter. The derived 
equations are essentially identical, but each approach has some advan- 
tages. Spinors make the the Petrov classification trivial once it has 
been shown that a tensor with the symmetries of the conformal tensor 
corresponds to a completely symmetric spinor, ^ abcd- The standard 
notation for the components of this tensor is 

#o = #oooo, #1 = #0001, ••• #4 = #1111. 
If £ is an arbitrary spinor then the equation 

*ABCD^ B C, C t D = 

is a homogeneous quartic equation with four complex roots, £i . The 
related real null vectors, Zf a = Cf Cf , are the four PNVs of Petrov. 
The spinor £ Q = Sq gives a PNV if \&o = 0. It is a repeated root and 
therefore it is the principal null vector of an AS spacetime if "Ji = as 
well. 

The Kerr (1963) letter presented the main results of my calculations 
but gave few detail^]. The methods that I used to solve the equations for 
AS spaces are essentially those used in "Exact Solutions" by Stephani 
et al. (2003), culminating in their equation (27.27). I will try to use the 
same notation as in that book since it is almost identical to the one I 
used in 1963, but I may get some of the signs wrong. Beware! 

f Robinson and Trautman also had a fairly natural complex tetrad approach. 

X I spent many years trying to write up this research but, unfortunately, I could 
never decide whether to use spinors or a complex tetrad, and thus it did not get 
written up until Kerr and Debney (1969). George Debney also collaborated with 
Alfred Schild and myself on the Kerr-Schild metrics in Debney et al. (1970). 
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Suppose that (e a ) = (m,fh, l,k) is a null tetrad, i.e., a set of four 
null vectors where the last two are real and the first two are complex 
conjugates. The corresponding dual forms are (w a ) = (fh,m, —k, —I) 
and the metric is 

ds 2 = 2(mm- kl) = 2{uj x uj 2 - w 3 w 4 )0 (1.3) 

The vector A; is a PNV and so its direction is uniquely defined, but 
the other directions are not. The form of the metric tensor in (|1.3[) is 
invariant under a combination of a null rotation (B) about k, a rotation 
(C) in the m Am plane and a Lorentz transformation (A) in the I A k 



plane. 












k' 


= k, 


m! = m + Bk, 


1' 


= l + Bfh + Bm + BBk, 


(1.4a) 


k' 


= k, 


m = e lC m, 


I' 


= 1, 


(1.4b) 


k' 


= Ak, 


ml — m. 


i 


= A-H. 


(1.4c) 



The most important connection form is 

r 4i = r 4 i a w a = m a k a .pdx 13 

The optical scalars of Ray Sachs for k are just the components of this 
form with respect to the u a . 



a = 


Fill 


= shear. 


P = 


r 4 i2 


= complex divergence 


K = 


r4i4 


= geodesy, 



The fourth component, r 4 i3, is not invariant under a null rotation 
about fc, 

r 4 i3 = r 4 i3 + Bp, 

and has no real geometric significance. It can be set to zero using an 
appropriate null rotation. Also, since k is geodesic and shearfree for AS 
metrics, both k and a are zero and therefore 

r 41 = ^ 2 . (i.5) 

If we use the simplest field equations, 

R44 = 2i? 4 l 4 2 = 0, i?41 = ^4112 — -^4134 = 0, Rn — 2i?4ii 3 = 0, 

f I personally hate the minus sign in this expression and did not use it in 1963, but 
it seems to have become standard. By the time I finish this article I am sure that 
I will wish I had stuck with all positive signs! 
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a total of 5 real equations, and the fact that the metric is AS, 

*0 = -2JZ4U1 = 0, 2*1 = -#4112 - #4134 = 0, 

then the most important of the second Cartan equations simplifies to 

dr 41 - r 4i a (ru + r 34 ) = R±i ab uj a a uo b = r^^ 2 a ^ 3 . (i.e) 

Taking the wedge product of (|1.6p with r 4 i and using (|1.5p . 

r 41 Adr 41 = o. (1.7) 

This was the key step in my study of these metrics but this result was 
not found in quite such a simple way. At first, I stumbled around using 
individual component equations rather than differential forms to look 
for a useful coordinate system. It was only after I had found this that 
I realized that using differential forms from the start would have short- 
circuited several days analysis. 

Equation (|1.7p is the integrability condition for the existence of com- 
plex functions, Q and II, such that 

r 41 = df/n, r 42 = dc/n. 

The two functions C, and its complex conjugate, £, will be used as (com- 
plex) coordinates. They are not quite unique since C can always be 
replaced by an arbitrary analytic function $ (Ca- 
using the transformations in (|1.4bj) and (ll.4c| . 

T vv = Ae tC T 41 = Ae lC d(/Il, => if = A^e'^Il. 

II' can be set to 1 by choosing Ae = II, and that is what I did in 1963, 
but it is also common to just use the C-transformation to convert it to 
a real function P, 

r 4i = pcu 2 = d(/p . (1.8) 

This is the derivation for two of the coordinates used in 1963. Since 
u)\k a = — > k(Q — 0, these functions, £, C, are constant along the 
PNV. 

The other two coordinates were very standard and were used by most 
people considering similar problems at that time. The simplest field 
equation is 

i? 44 = => kp = p\ 4 = p 2 , 
so that the real part of — p _1 is an affine parameter along the rays. This 
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was the obvious choice for the third coordinate, r, 

p- 1 = -(r + »S). 

There was no clear choice for the fourth coordinate, so u was chosen so 
that l a u a = 1, k a u a = 0, a pair of consistent equations. 
Given these four coordinates, the basis forms are 

uj 1 = m a dx a = -d(/Pp =(r- iE)d(/P, 
uo 2 = rh a dx a = -d(/Pp = (r + ffi)dC/P, 
u 3 = k a dx a = du + Ld( + LdQ, 
uj 4 = l a dx a = dr + Wd( + Wd( + Hlo 3 . 

where L is independent of R, and the coefficients E, W and H have still 
to be determined. 

When all this was substituted into the first Cartan equation, (|1 .42j) . 
and (|1.6|) . the simplest component of the second Cartan equation, (|1 .43jl . 
S and W were calculated as functions of L and its derivative^, 

2iY, = P 2 {dL~dL), d = d c -Ld u , 
W = -(r + iS)L„ + i9E. 

The remaining field equations, the "hard" ones, were more complicated, 
but still fairly straightforward to calculate. Two gave H as a function 
of a real "mass" function m(u, and certain functions of the higher 
derivatives of P and 10, 

, . mr + ME 

H=\K-r{\nP\ u - -g - - 2 . 

M = Y,K + P 2 Re[ddZ - 2L : „9E - Y.d u dL], 
K = 2P- 2 Re[d{d\nP - L, u )}, 

Finally, the first derivatives of the mass function, to, are given by the 
rest of the field equations, i?3i = and R33 = 0, 

d{m + iM) = 3(m + iM)L jU , (1.9a) 
B(m - iM) = 3(m - iM)L iU , (1.9b) 
[p- 3 {m + iM)], u = P[d + 2(91nP - L tU ]dI, (1.9c) 

f In Kerr (1963) f2, D and A were used instead of L, d and S, but the results were 
the same, mutatis mutandis. 

\ This expression for M was first published by Robinson et al. (1969). The corre- 
sponding expression in Kerr (1963) is for the gauge when P = 1. The same is true 
for equation jl.9t ). 



12 



R. P. Kerr 



where 

I = B{B\nP- Z „) + {B\nP- L :U f. (1.10) 

As was said in Kerr (1963), there are two natural choices that can 
be made to restrict the coordinates and simplify the final results. One 
is to rescale r so that P = 1 and L is complex, the other is to take 
L to be pure imaginary with P ^ 1. I chose to do the first since this 
gives the most concise form for M and the remaining field equations. It 
also gives the smallest group of permissible coordinate transformations, 
simplifying the task of finding all possible Killing vectors. The results 
for this gauge are 

M = lm(8ddL), (1.11a) 
3{m + iM) =3(m + iM)L tU , (1.11b) 
8{m-iM) = 3(m-iM)L„, (1.11c) 
d u [m - Re(dddL)} = \d u dL\ 2 . (l.lld) 

Since all derivatives of m are given, the commutators were calculated 
to see whether the field equations were completely integrable. This gives 
m as a function of the higher derivatives of L unless both S jU and L tUU 
are zero. As stated in Kerr (1963), if these are both zero then there 
is a coordinate system in which P and L are independent of u and 
m = cu + £), where c is a real constant. If it is zero then the metric 
is independent of u and therefore stationary. The field equations are 



V[V(lnP)] = c, V = P 2 d 2 /d(d(, (1.12a) 

M = 2SV(lnP)+VE, m = cu + A{(,(), (1.12b) 

cL = (A + iM) c , => VM = cE. (1.12c) 

We shall call these metrics quasi-stationary. 



In Kerr (1963) I stated that the solutions of these equations include 
the Kerr metric (for which c = 0). This is true but it is not how this 
solution was found. Furthermore, in spite of what many believe, its 
construction did not use the Kerr-Schild ansatz. 

1.3 Symmetries in Algebraically Special Spaces 

The field equations, (|1.9p or (jl.ll|) . are so complicated that some extra 
assumptions were needed to reduce them to a more manageable form. 
My next step in the hunt for physically interesting solutions was to 
assume that the metric is stationary. Fortunately, I had some tricks that 
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allowed me to find all possible Killing vectors without actually solving 
Killing's equation. 

The key observation is that any Killing vector generates a 1-parameter 
group which must be a subgroup of the group C of coordinate transfor- 
mations that preserve all imposed coordinate conditions. 

Suppose that {x* a , w*} is another set of coordinates and tetrad vectors 
that satisfy the conditions imposed in the previous section. If we restrict 
our coordinates to those that satisfy P = 1 then C is the group of 
transformations x — * x* for which 

C* = *(C), w 1 * = (l$d/*c)<A 

u* = |$ c |(« + 5(C,C), ^ = |$ ( |-V, 

r* = |$c|- 1 r J w 4 * = |$d<A 

and the basic metric functions, L* and m* , are given by 

L* = (\* C \/* C )[L-S C - i($cc/*c)(" + 5(C,C)], (l-13a) 
m ^ = |$ c |- 3 TO . (1.13b) 

Let <S be the identity component of the group of symmetries of a space. 
If we interpret these as coordinate transformations, rather than point 
transformations, then it is the set of transformations x — > x* for which 

9%p{x*) =g a p(x*). 

It can be shown that S is precisely the subgroup of C for which 

m*(x*) = m(x*), L*(x*) = L(x*)0 

Suppose that x — > x*(x,t) is a 1-paramcter group of motions, 

C* = ^(C;*), 

u* = \^\(u + T(CX;t), 
r* = 

Since x*(x; 0) = x, the initial values of ip and T are 

V(C;0) = C, T(c,C;0) = o. 

The corresponding infinitesimal transformation, K = K^d/dx 11 is 

~dx*»~ 



K» = 



dt 



t Note that this implies that all derivatives of these functions are also invariant, and 
so g a fs itself is invariant. 
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If we define 



«(0 



dip 



dt 



dT 



dt 



(1.14) 



then the infinitesimal transformation is 

K = ad^ + ad(- + Re(a^)(ud u — rd r ) + Vd u . 

Replacing $(£) with ip((;t) in (|1.13p . differentiating this w.r.t. t, and 
using the initial values for ip and T, K is a Killing vector if and only if 

V c + ±a cc r + KI+i(a ( - aj)L = 0, 
Krn + 3Re(a^)m = 0. 

The transformation rules for K under an element ($, S) of C are 

a* = $ c a, F* = \$ c \[V-Re(a c )S + KS}. 

Since a is itself analytic, if a ^ for a particular Killing vector then, 



if a 



then so is a*, and K is already simple without the (<£>, S) transformation 
being used. There are therefore two canonical types for K, 



Type 1 : K = Vd u , or Type 2 : K = d ( 



d, 



(1.15) 



These are asymptotically timelike and spacelike, respectively. 



1.4 Stationary solutions 

The obvious and easiest way to simplify the field equations was to as- 
sume that the metric was stationary. The Type 2 Killing vectors are 
asymptotically spacelike and so I assumed that (£ had a Type 1 Killing 

vector K = Vd u . The coordinates used in the last section assumed that 
l 

P = 1. If we transform to the more general coordinates where F ^ 1, 
using an A-transformation (jl.4cp with associated change in the (r, u) 
variables, we get 



k' = Ak, l' = A~Hy = A~ x r, 



K = Vd u = VAd u , = d u , if VA = 1. 

The metric can therefore be assumed independent of it, but P may not 
be constant. The basic functions, L,P and m are functions of (£,C) 



f Or any other analytic function of £ that one chooses 
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alone, and the metric simplifies to 

ds 2 = ds 2 + 2mr/(r 2 + E 2 )fc 2 , 
where the "base" metric, is 

{dsaf = 2{r 2 + -£ 2 )p- 2 d(d( - 2l„k, 



lo = dr + i(E x d( -T, rf() 



ME 



(r 2 + E 2 ) 



(1.16) 

(1.17a) 
(1.17b) 



Although the base metric is flat for Schwarzschild it is not so in gen- 
eral. E, K and M are all functions of the derivatives of L and P. 



K = 2V 2 lnP, 

72 rtli 



(1.18) 



E = P 2 lm(L^), 

M = Y,K + V 2 E, V 2 = P 2 d c d^, 

The mass function, m,and M are conjugate harmonic functions, 

m c = -zM c , m^ = +iM^, (1.19) 

and the remaining field equations are 

V 2 K = 0, V 2 Af = 0. (1.20) 

If to is a particular solution of these equations then so is to+too where too 
is an arbitrary constant. The most general situation where the metric 
splits in this way is when P, L and M arc all independent of u but 
to = cu + A((, £). The field equations for these are given in p. 120 and 
in Kerr (1963). We can state this as a theorem: 

Theorem 2 If ds% is any stationary (diverging) algebraically special 
metric, or more generally a solution of i LIS]) , then so is 



ds n 



2mor 
r 2 + E : 



k 2 , 



where too is an arbitrary constant. These are the most general diverging 
algebraically special spaces that split in this way. 

These are "generalized Kerr-Schild" metrics with base spaces ds 2 , that 
are not necessarily flat. 

The field equations for stationary AS metrics are certainly simpler 
than the original ones, (|1.9|) . but they are still PDEs, not ODEs, and 
their complete solution is unknown. 
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1.5 Axial symmetry 

We are getting close to Kerr. At this point I assumed that the metric 
was axially symmetric as well as stationary. I should have revisited the 
Killing equations to look for any Killing vector (KV) that commutes 
with the stationary one, d u . However, I knew that it could not also be 
Type 10 and therefore it must be Type 2. It seemed fairly clear that 
it could be transformed to the canonical form i(d^ — d^) (= d y where 
£ = x + iy) or equivalently i(C@C ~ C^f ) (— dtf> m polar coordinates where 
( = Re 1 ^). I was getting quite eager at this point so I decided to just 
assume such a KV and see what turned urjj 

From the first equation in ()1 .20|) . the curvature V 2 (lnP) of the 2- 
metric P _2 d£d£ is a harmonic function, 

V 2 lnP = P 2 (lnP), cf = F(0 + F(C), 

where F is analytic. There is only one known solution of this equation 
for F not a constant, 

P 2 = Po(C + C) 3 , V 2 lnP = -^P (C + C), (1-21) 

where Pq is an arbitrary constant. The mass function m is then constant 
and the last field equation, V 2 A/ = 0, can be solved for L. The final 
metric is given in Kerr and Debney (1969), equation (6.14), but it is not 
worth writing out here since it is not asymptotically flat. 

If (£ is to be the metric for a localized physical source then the null 
congruence should be asymptotically the same as Schwarzschild. P(C) 
must be regular everywhere, including at infinity, and must therefore be 
constant, 

\K = PP CC - - P C P C - = R = ±P 2 , (say). (1.22) 

As was shown in Kerr and Debney (1969), the appropriate Killing 
equations for a K that commutes with K are 

2 1 

K = adQ + ad^, a = a(C), 
KL = -a c L, KE = 0, (1.23) 
KP = Re(a c )P, Km = 0. 

f If it were it would be parallel to the stationary KV and therefore a constant 
multiple of it. 

| All possible symmetry groups were found for diverging AS spaces in George C. 
Debney's Ph.D. thesis. My 1963 expectations were confirmed there. 
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I do not remember the choice I made for the canonical form for K in 

2 

1963, but it was probably d y . The choice in Kerr and Debney (1969) 
was 

a = K, K = ^ - 

and that will be assumed here. For any function f(C,C)i 

K/ = => f(£,0 = g(Z), where Z = & 

Now Re(a ) f) = 0, and therefore 

KP = 0, => P = P(Z), 

2 

and therefore K is given by 

\K = P 2 (In P) <a = PP a - P c P f = Z P = Z + Z , 

after a <i>(C)-coordinate transformation. Note that the form of the metric 
is invariant under the transformation 

r = A Q r*, u = A7 1 u*, C = A (* , 

_ 2 o (1-24) 

Z Q = A Zg, m = A m , 

where Aq is a constant, and therefore Zq is a disposable constant. We 
will choose it later. 

The general solution of (ll.23|) for L and S is 

L = iCP~ 2 B(Z), Z = ZB' - (1- ZaP-^B, 

where B' = dB/dZ. The complex "mass", m + iM, is an analytic 
function of ( from (|1.19[) . and is also a function of Z from (|1.23[) . It 
must therefore be a constant, 

m + iM = (J® = mo + iAfp. 
Substituting this into (| 1 . 1 8(1 . the equation for E, 
Eif + V 2 E = M = M — ► 
P 2 [ZYI' + E'] + 2Z E = M Q . 
The complete solution to this is 

Z = C + ?—^l-a + C 2 \nZ}, 

Zj + Zq 

where {Co, a, C 2 } are arbitrary constants. This gave a four-parameter 
metric when these known functions are substituted into (|1.16p and (|1.17p . 
However, if C 2 is nonzero then the final metric is singular at r = 0. It 
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was therefore omitted in Kerr (1963). The "imaginary mass" is then 
M = 2ZqCq and so Cq is a multiple of the NUT parameter. It was 
known in 1963 that the metric cannot be asymptotically flat if this is 
nonzero and so it was also omitted. The only constants retained were 
mo, a and Z . When a is zero and Z is positive the metric is that of 
Schwarzschild. It was not clear that the metric would be physically in- 
teresting when a 0, but if it had not been so then this whole exercise 
would have been futile. 

The curvature of the 2-metric, 2P~~ 2 d(d(, had to have the same sign 
as Schwarzschild if the metric was to be asymptotically fiat, and so 
Zq = +Pq- The basic functions in the metric then became 

Z Q = Pi, P = (C + Pi, m = m , M = 0, 

r • CC — Zq 

CC + Zq 

The metric was originally published in spherical polar coordinates. 
The relationship between these and the (C, C) coordinates is 

C = Pocotf e**. 

At this point we choose Aq in the transformation (|1.24|) so that 

2Pi = 1, =>■ k = du + a sin 2 0d(j) 
Recalling the split of (fTTTC|) and ifEIT]) . 

ds 2 = ds 2 + 2mr/(r 2 + a 2 cos 2 6)k 2 (1.25) 

where m = mo, a constant, and 

dsl =(r 2 + a 2 cos 2 9){d9 2 + sin 2 '6d<j) 2 ) 

— (2dr + du — a sin Odtfi) (du + a sin " ' 



This is the original form of Kerr (1963), except that u has been replaced 
by — u to agree with current conventions, and a has been replaced with 
its negativqjj. 

Having found this fairly simple metric, I was desperate to see whether 
it was rotating. Fortunately, I knew that the curvature of the base 
metric, ds^, was zero, and so it was only necessary to find coordinates 
where this was manifestly Minkowskian. These were 

(r + ia)e l ^sm9 — x + iy, r cos# = z, r + u = —t. 



\ We will see why later. 
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This gives the Kerr-Schild form of the metric, 

ds 2 = dx 2 + dy 2 + dz 2 - dt 2 + —. ^[dt + -dz 

r 4 + a z r ^ 27) 

+ -5-. — o(xdx + ydy) — — ^(xdy - ydx)] 2 . 

r z _|_ a z r z _|_ a z 

where the surfaces of constant r are confocal ellipsoids of revolution 
about the z-axis, 

r z + a A r A 

Asymptotically, of course, r is just the distance from the origin in the 
Minkowskian coordinates, and the metric is clearly asymptotically flat. 



Angular momentum 

After the metric had been put into its Kerr-Schild form I went to Alfred 
Schild and told him I was about to calculate the angular momentum of 
the central body. He was just as excited as I was and so he joined me in 
my office while I computed. We were both heavy smokers at that time, 
so you can imagine what the atmosphere was like, Alfred puffing away 
at his pipe in an old arm chair, and myself chain-smoking cigarettes at 
my desk. 

The Kerr-Schild form is particularly suitable for calculating the phys- 
ical parameters of the solution. My PhD thesis at Cambridge was enti- 
tled "Equations of Motion in General Relativity" . It had been claimed 
previously in the literature that it was only necessary to satisfy the mo- 
mentum equations for singular particles to be able to integrate the EIH 
quasi-static approximation equations at each order. One thing shown in 
my thesis was the physically obvious fact that the angular momentum 
equations were equally important. Some of this was published in Kerr 
(1958) and (1960). Because of this previous work I was well aware how 
to calculate the angular momentum in this new metric. 

It was first expanded in powers of J? -1 , where R = x 2 + y 2 + z 2 is the 
usual Euclidean distance, 



ds 2 =dx 2 + dy 2 + dz 2 - dt 2 + ^(dt + dR) 2 

R 

-(xdy-ydx)(dt + dR) + 0{R- 3 ) 



(1.29) 



R 3 

Now, if x^ — ► x^+a^ is an infinitesimal coordinate transformation, then 
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ds 2 — > (is 2 + 2da^dx^ . If we choose 



a^efe^ = -^(xdy - ydx) => 



Ida^dx 11 = —4m — — (xdy — ydx)dR, 



4am 

then the approximation in (| 1 . 29[) simplifies to 



ds 2 =dx 2 + dy 2 + dz 2 - dt 2 + —(dt + dR) 2 

R (1.30) 



4ma 



(xdy - ydx)dt + 0(R~ 6 ) 



R 3 

The leading terms in the linear approximation for the gravitational field 
around a rotating body were well known (for instance, see Papapetrou 
(1974) or Kerr (I960)). The contribution from the angular momentum 
vector, J, is 

4R- z e ljk J l x°dx k dt. 

A comparison of the last two equations showed that the physical param- 
eters were 

Mass = m, J = (0, 0, ma)[j] 

When I turned to Alfred Schild, who was still sitting in the arm-chair 
smoking away, and said "Its rotating!" he was even more excited than I 
was. I do not remember how we celebrated, but celebrate we did! 

Robert Boyer subsequently calculated the angular momentum by com- 
paring the known Lense-Thirring results for frame dragging around a 
rotating object in linearized relativity with the frame dragging for a cir- 
cular orbit in a Kerr metric. This was a very obtuse way of calculating 
the angular momentum since the approximation (|1.30|) was the basis for 
the calculations by Lense and Thirring, but it did show that the sign 
was wrong in the original paper! 



1.6 Singularities and Topology 

The first Texas Symposium on Relativistic Astrophysics was held in 
Dallas December 16-18, 1963, just a few months after the discovery of the 
rotating solution. It was organized by a combined group of Relativists 

f Unfortunately, I was rather hurried when performing this calculation and got the 
sign wrong. This is why the sign of the parameter a in Kerr(1963) is different to 
that in all other publications, including this one. This way of calculating J was 
explained at the First Texas Symposium at the end of 1963, see Kerr(1965), but 
I did not check the sign at that time. 
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and Astrophysicists and its purpose was to try to find an explanation for 
the newly discovered quasars. The source 3C273B had been observed 
in March and was thought to be about a million million times brighter 
than the sun. 

It had been long known that a spherically symmetric body could col- 
lapse inside an event horizon to become what was to be later called a 
black hole by John Wheeler. However, the Schwarzschild solution was 
non-rotating and it was not known what would happen if rotation was 
present. I presented a paper called "Gravitational collapse and rota- 
tion" in which I outlined the Kerr solution and said that the topological 
and physical properties of the event horizon may change radically when 
rotation is taken into account. It was not known at that time that Kerr 
was the only possible stationary solution for such a rotating black hole 
and so I discussed it as an example of such an object. 

Although this was not pointed out in the original letter, Kerr (1963), 
the geometry of Kerr is even more complicated than the Kruskal exten- 
sion of Schwarzschild. The metric is everywhere nonsingular, except on 
the ring 

2 = 0, x 2 + y 2 = a 2 . 

The Weyl scalar, R a bcd,R abcd ^ oo near these points and so they are true 
singularities, not just coordinate ones. Furthermore, this ring behaves 
like a branch point in the complex plane. If one travels on a closed curve 
that threads the ring the initial and final metrics are different: r changes 
sign. Equation (|1.28p has one nonnegative root for r 2 , and therefore two 
real roots, r±, for r. These coincide where r 2 — 0, i.e., on the disc D 
bounded by the ring singularity 

D: z = 0, x 2 +y 2 <a 2 . 

The disc can be taken as a branch cut for the analytic function r. We 
have to take two spaces, E\ and E% with the topology of R A less the disc 
D. The points above D in E\ are joined to the points below D in Ei 
and vice versa. In E\, r > and the mass is positive at infinity; in E^^ 
r < and the mass is negative. The metric is then everywhere analytic 
except on the ring. 

It was trivially obvious to everyone that if the parameter a is very 
much less than m then the Schwarzschild event horizon at r = 2m will 
be modified slightly but cannot disappear. For instance, the light cones 
at r — m in Kerr all point inwards for small a. Before I went to the 
meeting I had calculated the behaviour of the time like geodesies up and 
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down the axis of rotation and found that horizons occurred at the points 
on the axis in E\ where 

r 2 — 2mr + a 2 = 0, r = \z\. 

but that there are no horizons in E2 where the mass is negative. In 
effect, the ring singularity is "naked" in that sheet. 

I made a rather hurried calculation of the two event horizons in E\ 
before I went to the Dallas Symposium and claimed incorrectly there, 
Kerr (1964), that the equations for them were the two roots of 

r 4 - 2mr 3 + a 2 z 2 = 0, 

whereas z 2 should be replaced by r 2 in this and the true equation is 

r 2 - 2mr + a 2 = 0. 

I attempted to calculate this using inappropriate coordinates and as- 
suming that the equation would be: "ip{r, z) is null for some function of 
both r and z" . I did not realize that this function depended only on r. 

The Kerr-Schild coordinates arc a generalization of the Eddington 
Finkclstcin coordinates for Schwarzschild. For the latter future-pointing 
radial geodesies are well behaved but not those traveling to the past. 
Kruskal coordinates were designed to handle both. Similarly for Kerr, 
the coordinates given here only handle ingoing curves. This metric is 
known to be Type D and therefore it has another set of Debever-Penrose 
vectors and an associated coordinate system for which the outgoing 
geodesies are well behaved, but not the ingoing ones. 

This metric consists of three blocks, outside the outer event hori- 
zon, between the two horizons and within the inner horizon (at least 
for m < a, which is probably true for all existing black holes). Just 
as Kruskal extends Schwarzschild by adding extra blocks, Boyer and 
Lindquist (1967) and Carter (1968) independently showed that the max- 
imal extension of Kerr has a similar proliferation of blocks. However, 
the Kruskal extension has no application to a real black hole formed by 
the collapse of a spherically symmetric body and the same is true for 
Kerr. In fact, even what I call E 2 , the sheet where the mass is negative, 
is probably irrelevant for the final state of a collapsing rotating object. 

Ever since this metric was first discovered people have tried to fit an 
interior solution. One morning during the summer of 1964 Ray Sachs 
and myself decided that we would try to do so. Since the original form is 
useless and the Kerr-Schild form was clearly inappropriate we started by 
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transforming to the canonical coordinates for stationary axisymmctric 
solutions. 

In Papapetrou (1966) there is a very elegant treatment of stationary 
axisymmetric Einstein spaces. He shows that if there is a real non- 
singular axis of rotation then the coordinates can be chosen so that 
there is only one off-diagonal component of the metric. We call such a 
metric quasi-diagonalize. All cross terms between {dr, dO} and {dt, d(j>} 
can be eliminated by transformations of the type 

dt' =dt + Adr + Bd6, d(j>' = d<f) + Cdr + DdO. 

where the coefficients can be found algebraically. Papapetrou proved 
that dt' and d(f>' are perfect differentials if the axis is regulaJJ 

Ray and I calculated the coefficients A. . . D and transformed the met- 
ric to the Boyer-Lindquist form, 

2mr 

dt — > dt H t — dr 

A 

d<f> — > —d<fi + -rdr 

A = r 2 — 2mr + a 2 , 
S = r 2 + a 2 cos 2 6, 

where, as before, u = —(t+r). The right hand sides of the first two equa- 
tions are clearly perfect differentials as the Papapetrou analysis showed. 
The full Boyer-Lindquist form of the metric is 

£ A 
ds 2 =~rdr 2 [dt - asm 2 6d(j)} 2 + 

2fl (1-31) 
^ + h .^-[(r 2 +a 2 )d^-adt] 2 , 

after some tedious analysis that used to be easy but now seems to require 
an algebraic package such as Maple. 

Having derived this canonical form, we studied the metric for at least 
ten minutes and then decided that we had no idea how to introduce a 
reasonable source into a metric of this form, and probably would never 
have. Presumably those who have tried to solve this problem in the 
last 43 years have had similar reactions. Soon after this failed attempt 
Robert Boyer came to Austin. He said to me that he had found a new 
quasi-diagonalized form of the metric. I said "Yes. It is the one with 

f It is shown in Kerr and Weir (1976) that if the metric is also algebraically special 
then it is quasi-diagonalize precisely when it is Type D. These metrics are the 
NUT parameter generalization of Kerr. 
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the polynomial r 2 — 2mr + a 2 " but for some reason he refused to believe 
that we had also found this form. Since it did not seem a "big deal" 
at that time I did not pursue the matter further, but our relations were 
hardly cordial after that. 

One of the main advantages of this form is that the event horizons 
can be easily calculated since the inverse metric is simple. If f(r, 9) = 
is a null surface then 

A(r)/ >rr + lee = 0, 

and therefore A < 0. The two event horizons are the surfaces r = r± 
where the parameters r± are the roots of A = 0, 

A = r 2 — 2mr + a 2 = (r — r + )(r — r_). 

If a < to there are two distinct horizons between which all time-like 
lines point inwards; if a = to there is only one event horizon; and for 
larger a the singularity is bare! Presumably, any collapsing star can only 
form a black hole if the angular momentum is small enough: a < to. This 
seems to be saying that the body cannot rotate faster than light, if the 
final picture is that the mass is located on the ring radius a. However, 
it should be remembered that this radius is purely a coordinate radius, 
and that there is no way that the final stage of such a collapse is that 
all the mass is located at the singularity. 

The reason for the last statement is that if the mass where to end on 
the ring then there would be no way to avoid the second asymptotically 
flat sheet where the mass appears negative. I do not believe that the 
star opens up like this along the axis of rotation. If we remember that 
the metric is discontinuous across the disc bounded by the singular ring 
then it is quite possible that a well-behaved finite source could be put 
between z = 0±, \R\ < cp. This would correspond to the surface of the 
final body being r = in Boyer-Lindquist coordinates, say, but where 
the interior corresponds to r < 0.. The actual surface may be more 
complicated than this but I am quite sure that this is closer to the final 
situation than that the matter all collapses onto the ring. 

What I believe to be more likely is that the inner event horizon never 
actually forms. As the body continues to collapse inside its event horizon 
it spins faster and faster so that the geometry in the region between its 
outer surface and the outer event horizon approaches that between the 
two event horizons for Kerr. The surface of the body surface will appear 

f This has been done using S— functions, but I am thinking more of a nonsingular 
source where the distance between the two sides of the disk is nonzero. 
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to be asymptotically null. The full metric may not be geodesically com- 
plete. Many theorems have been claimed stating that a singularity must 
exist if certain conditions arc satisfied, but they all make assumptions 
that may not be true for collapse to a black hole. Furthermore, these as- 
sumptions are often (usually?) unstated or unrecognised, and the proofs 
are dependent on other claims/theorems that may not be correct. 

These are only two of a very large range of possibilities for the interior. 
What happens after the outer horizon forms is still a mystery after 
more than four decades. It is also the main reason why I said at the 
end of Kerr (1963) that "It would be desirable to calculate an interior 
solution. . . " . This statement has been taken by some to mean that I 
thought the metric only represented a real rotating star. This is untrue 
and is an insult to all those relativists of that era who had been looking 
for such a metric to see whether the event horizon of Schwarzschild would 
generalise to rotating singularities. 

The metric was known to be Type D with two distinct geodesic and 
shearfree congruences from the moment it was discovered. This means 
that if the other is used instead of k then the metric must have the 
same form, i.e., it is invariant under a finite transformation that reverses 
"time" and possibly the axis of rotation in the appropriate coordinates. 
This also meant that there is an extension that is similar to the Kruskal- 
Szekeres extension of Schwarzschild. Both Boyer and Lindquist (1967) 
and a fellow Australasian, Brandon Carter (1968), solved the problem of 
constructing the maximal extensions of Kerr, and even that for charged 
Kerr. These are mathematically fascinating and the latter paper is a 
beautiful analysis of the problem, but the final result is of limited phys- 
ical significance. 

Brandon Carter's (1968) paper was one of the most significant papers 
on the Kerr metric during the mid-sixties for another reason. He showed 
that there is an extra invariant for geodesic motion which is quadratic 
in the momentum components: J = X a t,v a v b where X a b is a Killing 
tensor, A( afc;c ) = 0. This gave a total of four invariants with the two 
Killing vector invariants and |v| 2 itself, enough to generate a complete 
first integral of the geodesic equations. 

Another significant development was the "proof" that this is the only 
stationary metric with a simply connected bounded event horizon, i.e., 
the only possible black hole. Many contributed to this, including Steven 
Hawking (1972), Brandon Carter (1971) and another New Zealander, 
David Robinson (1975). The subsequent work in this area is discussed 
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by David in an excellent article in this book and so I will not pursue 
this any further here. 



1.7 Kerr Schild metrics 

One morning during the fall semester 1963, sometime before the First 
Texas Symposium, I tried generalizing the way that the field equations 
split for the Kerr metric by setting 

ds 2 — dsl + -= -rrrrfc 2 . 

The base metric ds^ was to be an algebraically special metric with m — 0. 
From an initial rough calculation this had to be flat. Also, it seemed 
that the coordinates could be manipulated so that dL = Lq — LL U = 
and that the final metric depended on an arbitrary analytic function of 
the complex variable £. At this point I lost interest since the metric had 
to be singular at the poles of the analytic function unless this function 
was quadratic and the metric was therefore Kerr. 

Sometime after the Texas Symposium Jerzy Plebanski visited Austin. 
Alfred Schild gave one of his excellent parties for Jerzy during which I 
heard them talking about solutions of the Kerr-Schild typc0, 
dsQ + hk 2 , where the first term is flat and k is any null vector. I com- 
mented that there might be some algebraically special spaces with this 
structure depending on an arbitrary function of a complex variable but 
that this had not been checked. 

At this point Alfred and I retired to his home office and calculated the 
simplest field equation, R a bk a k b — 0. To our surprise this showed that 
the null vector had to be geodesic. We then calculated k[ a R b ] pq [ c k d ]k p k q , 
found it to be zero and deduced that all metrics of this type had to be 
algebraically special. This meant that all such spaces with a diverging 
congruence might already be known. We checked my original calcula- 
tions next day and found them to be correct. 

As was stated in Theorem [2j m is a unique function of P and L unless 
there is a canonical coordinate system where m is linear in u and {L, P} 
are functions of {£, alone. If the base space is flat then m :U = c = 
and the metric is stationary. The way these Kerr-Schild metrics were 
found originally was by showing that in a coordinate system where P = 1 
the canonical coordinates could be chosen so that dL = 0. Transforming 



f This name came later, of course. 
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from these coordinates to ones where P^l and d u is a Killing vector, 

P cc = 0, L = P~ 2 <HC), (1.32) 

where 0(C) is analytic. From the first of these P is a real bilinear function 
of ( and therefore of 

P = KC + q( + ?C + c 

This can be simplified to one of three canonical forms, P = 1, 1 ± £C by 
a linear transformation on £. We will assume henceforth that 

p = l + cc 

The only problem was that this analysis depended on results for alge- 
braically special metrics and these had not been published and would 
not be for several years. We had to derive the same results by a more 
direct method. The metric was written as 

ds 2 = dx 2 + dy 2 + dz 2 - dt 2 + hk 2 , (1.33a) 
k = (du + Yd( + Yd( + YYdv)/{\ + YY), (1.33b) 

where Y is the old coordinate £ used in (|1.32|fj l and 

u = z + t, v — t — z, C = x + ijjj]- 

The tetrad used to calculate the field equations was defined naturally 
from the identity 

ds 2 = (d( + Ydv)(d( + Ydv) - (dv + hk)k 

Each of these spaces has a symmetry which is also a translational 
symmetry for the base Minkowski space, ds 2 ). The most interesting 
situation is when this is time-like and so it will be assumed that the 
metric is independent of t = |(it + v). 

If (f>(Y) is the same analytic function as in (|1.32|) then Y is determined 
as a function of the coordinates by 

Y 2 C + 2zY - C + <j){Y) = 0. (1.34) 

and the coefficient of k 2 in p. 33ft is 

h = 2mRe(2Y c ), (1.35) 

f Y is the ratio of the two components of the spinor corresponding to fc. 
| Note that there certain factors of \/2 have been omitted to simplify the results. 
This does lead to a factor 2 appearing in 111,350 , 
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where m is a real constant. Differentiating (|1.34[) with respect to £ gives 
Y c = (2Y( + 2z + <f>')- 1 . (1.36) 
The Weyl spinor invariant is given by 

* 2 = c mY" c 3 , 

where cq is some power of 2, and the metric is therefore singular precisely 
where Y is a repeated root of its defining equation, (|1.34|) . 

If the k-lines are projected onto the Euclidean 3-space t = with 
{x, y, z) as coordinates so that ds 2 E = dx 2 + dy 2 + dz 2 , then the perpen- 
dicular from the origin meets the projected k-line at the point 

d>-Y 2 4> Ycf> + Y4> 

and the distance of the line from the origin is 

101 



D = 



i + yy' 



a remarkably simple result. This was used by Kerr and Wilson (1978) to 
prove that unless <fi is quadratic the singularities are unbounded and the 
spaces are not asymptotically flat. The reason why I did not initially 
take the general Kerr-Schild metric seriously was that this was what I 
expected. 

Another point that is easily calculated is Zq where the line meets the 
plane z = 0, 

4> + y 2 4> 

Z ° : C= 1-(YY) 2 ' Z = °' 
The original metric of this type is Kerr where 

If 4>(Y) is any other quadratic function then it can be transformed to the 
same value by using an appropriate Euclidean rotation and translation 
about the t-axis. The points Fq and Zq are the same for Kerr, so that 
Fq lies in the z-plane and the line cuts this plane at a point inside the 
singular ring provided \ Y\ ^ 1. The lines where \Y\ = 1 are the tangents 
to the singular ring lying entirely in the plane z = outside the ring. 
When a — > the metric becomes Schwarzschild and all the F-lines pass 
through the origin. 
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When 4>(Y) = -2iaY (fl~34|) becomes 

Y 2 C + 2(z-ia)Y-C = 0. 
There are two roots, Y± and Y% of this equation, 

rC r 

y x = 7 — ■ ,; — — , 2y u = 



(z + r)(r — ia) ' r 4 + a 2 z 2 

_ < _ r 3 + zarz 

-12 — 7 T7 : — r, ^-12 C — i t — 9 

(z — r)(r + za) + a z 2 z 

where r is a real root of (| 1 . 28[) . This is a quadratic equation for r 2 with 
only one nonnegative root and therefore two real roots differing only 
by sign, ±r. When these are interchanged, r <-> — r, the corresponding 
values for Y are also swapped, Y\ *r+Yz. 

When Y2 is substituted into the metric then the same solution is re- 
turned except that the mass has changed sign. This is the other sheet 
where r has become negative. It is usually assumed that Y is the first 
of these roots, Y\. The coefficient h of k 2 in the metric, (|1.33|) . is then 



h = 2mRe(2F c ) 



2mr A 



. * + a 2 z 2 

This gives the metric in its KS form. (|I.27p . 

The results were published in two places, Kerr and Schild ( 1965a, b). 
The first of these was a talk that Alfred gave at the Galileo Centennial 
in Italy, the second was an invited talk that I gave, but Alfred wrote, at 
the Symposium on Applied Mathematics of the American Mathematical 
Society, April 25 1964. The manuscript had to be presented before the 
conference so that the participants had some chance of understanding 
results from distant fields. We stated on page 205 that 

"Together with their graduate student, Mr. George Debney, the authors have 
examined solutions of the nonvacuum Einstein-Maxwell equations where the 
metric has the form (2.ljJ Most of the results mentioned above apply to this 
more general case. This work is continuing.". 



1.8 Charged Kerr 

What was this quote referring to? When we had finished with the Kerr- 
Schild metrics, we looked at the same problem with a nonzero electro- 
magnetic field. The first stumbling block was that Rabk a k b = no 

f Equation (1 1 . 33 ^ in this paper. It refers to the usual Kerr-Schild ansatz. 
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longer implied that the k-lines are geodesic. The equations were quite 
intractable without this and so it had to be added as an additional as- 
sumption. It then followed that the principal null vectors were shearfree, 
so that the metrics had to be algebraically special. The general forms 
of the gravitational and electromagnetic fields were calculated from the 
easier field equations. The E-M field proved to depend on two functions 
called A and 7 in Debney, Kerr and Schild (1969). 

When 7 = the "difficult" equations are linear and similar to those for 
the purely gravitational case. They were readily solved giving a charged 
generalization of the original Kerr-Schild metrics. The congruences arc 
the same as for the uncharged metrics, but the coefficient of k 2 is 

/i = 2mRe(2Y; c )-|i/f |2F C | 2 . (1.37) 

where ip(Y) is an extra analytic function generating the electromagnetic 
field. This is best expressed through a potential, 

/ = \F llv dx ll dx v = -da, 

a = -P(tpZ + $Z)k - \{ X dY + xdY), 

where 



/ 



p- 2 ip(Y)dY, 



Y being kept constant in this integration. 

The most important member of this class is charged Kerr. For this, 

r 4 + a 2 z 2 

Asymptotically, r = R, k = dt — dR, a radial null-vector and Y = 
tan(i#)e*^. If the analytic function ip(Y) is nonconstant then it must be 
singular somewhere on the unit sphere and so the gravitational and elec- 
tromagnetic fields will be also. The only physically significant charged 
Kerr-Schild is therefore when ip is a complex constant, e + ib. The imag- 
inary part, b, can be ignored as it gives a magnetic monopole, and so we 
are left with tjj = e, the electric charge, 

2777.7*^ ■ G^V^ Z 

ds 2 = dx 2 + dy 2 + dz 2 - dt 2 H -. [dt + -dz 

r 4 + a z z z r 

+ 2 7 2 (xdx + ydy) - ° (xdx - ydy)] 2 

The electromagnetic potential is 

- ; a(xdy — ydx) 



(1.39) 



er 

a = 



r 4 + a 2 z 2 



dt - 



r z + a z 
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where a pure gradient has been dropped. The electromagnetic field is 

(Fxt iFy Z7 Fyt iF ZXl F z t iF X y) 

Q 

(r 2 + iaz) A 

In the asymptotic region this field reduces to an electric field, 
and a magnetic field, 

H = ^(3x Zl 3yz,3z 2 - R 2 ). 

This is the electromagnetic field of a body with charge e and magnetic 
moment (0,0, ea). The gyromagnetic ratio is therefore ma/ea — m/e, 
the same as that for the Dirac electron. This was first noticed by Bran- 
don Carter and was something that fascinated Alfred Schild. 

This was the stage we had got to before March 1964. We were unable 
to solve the equations where the function 7 was nonzero so we enlisted 
the help of our graduate student, George Debney. Eventually we realized 
that we were unable to solve the more general equations and so we 
suggested to George that he drop this investigation. He then tackled 
the problem of finding all possible groups of symmetries in diverging 
algebraically special spaces. He succeeded very well with this, solving 
many of the ensuing field equations for the associated metrics. This 
work formed the basis for his PhD thesis and was eventually published 
in Kerr and Debney (1970). 

In Janis and Newman (1965) and Newman and Janis (1965) the au- 
thors defined and calculated multipole moments for the Kerr metric, 
using the Kerr-Schild coordinates as given in Kerr (1963). They then 
claimed that this metric is that of a ring of mass rotating about its 
axis of symmetry. Unfortunately, this cannot be so because the metric 
is multivalued on its symmetry axis and is consequently discontinuous 
there. The only way that this can be avoided is by assuming that the 
space contains matter on the axis near the centre. As was acknowledged 
in a footnote to the second paper, this was pointed out to the authors 
by the referee and myself before the paper was published, but they still 
persisted with their claim. 
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1.9 Newman's construction of the Kerr Newman metric 

Newman knew that the Schwarzschild, Reissner-Nordstrom and Kerr 
metrics all have the same simple form in Eddington-Finkelstein or Kerr- 
Schild coordinates (see eq. I1.33[) . Schwarzschild and its charged general- 
isation have the same null congruence, k — dr — dt; only the coefficient 
h is different, 



For these the complex divergence of the underlying null congruence is 



Newman hoped to find a charged metric with the same congruence as 
Kerr but with h generalised to something like the Reissner-Nordstrom 
form with e 2 /r 2 replaced by e 2 p 2 . This does not quite work since p is 
complex for Kerr so he had to replace p 2 with something real. 

There are many real rational functions of p and p that reduce to p 2 
when p is real, so he wrote down several possibilities and distributed 
them to his graduate students. Each was checked to see whether it was 
a solution of the Einstein-Maxwell equations. The simplest, p 2 — ■> pp, 
worked! The appropriate electromagnetic field was then calculated, a 
non-trivial problem. 

The reason that this approach was successful has nothing to do with 
"complexifying the Schwarzschild and Reissner-Nordstrom metrics" by 
some complex coordinate transformation, as stated in the original pa- 
pers. It works because all these metrics are of Kerr-Schild form and the 
general Kerr-Schild metric can be charged by replacing the uncharged h 
with its appropriate charged version, h = 2mRe(2y^) — |?/;| 2 |2Y^| 2 — ► 
2mRe(p) — e 2 pp, without changing the congruence. 

The charged solution was given in Newman et al. (1965). They claimed 
that the metric can be generated by a classical charged rotating ring. As 
in the previous paper Newman and Janis (1965), it was then admitted 
in a footnote that the reason why this cannot be true had already been 
explained to them. 

1.10 Appendix: Standard Notation 

Let {e a } and {uj a } be dual bases for tangent vectors and linear 1-forms, 
respectively, i.e., us a (et,) = 6%. Also let g a b be the components of the 




2m 




P = P = 1/r. 
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metric tensor, 

ds 2 = g ab u> a u> a , g a b = e a ■ e b . 

The components of the connection in this frame are the Ricci rotation 
coefficients, 

pa a fi v p ps 

1 be — ^ ii\veb e c , L a bc — gas 1 - bci 

The commutator coefficients D a b c = —D a cb are defined by 

[e b ,e c ] = D a bc e a , where [u, v](f) = u(v(f)) - v(u(f)). 
or equivalently by 

dcu a =DVAw c . (1.40) 

Since the connection is symmetric, D a b c = — ^ a [bc], an d since it is 
metrical 

Tabc = \{g a b\c + 9ac\b — 9bc\a + -Dfcac + Dcab ~ D a bc), 
abc — yam 1 be: U abc — yam*-J be- 

If it is assumed that the g a b are constant, then the connection compo- 
nents are determined solely by the commutator coefficients and therefore 
by the exterior derivatives of the tetrad vectors, 

I" abc ~^[J-)bac ~l~ T) ca b Dabc)- 

The components of the curvature tensor are 

Q a bcd = r a M | C — r a bc | d + T e bdF a ec — r e b c r a e( j — D e C dV a b e - (i-4i) 

We must distinguish between the expressions on the right, the @ a bcd, 
and the curvature components, R a bcd, which the N-P formalism treat 
as extra variables, their (^>i). 

A crucial factor in the discovery of the spinning black hole solutions 
was the use of differential forms and the Cartan equations. The connec- 
tion 1-forms F a b are defined as 

ra j^a , c 

6—1 bc^ ■ 

These are skew-symmetric when g a b\ c = 0, 

Tfra F a fr, r a fr — (/ ac T a - 

The first Cartan equation follows from (|1.40p . 

duj a + T a b uo b = 0. (1.42) 
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The curvature 2-forms are defined from the second Cartan equations, 
&\ = dT\ + T a c A r c 6 = \R a bcd u c u d . (1.43) 
The exterior derivative of (|1.42j) gives 

©\ac/ = o e a [bcd] = o, 

which is just the triple identity for the Riemann tensor, 

R a [bcd] = 0. (1.44) 
Similarly, from the exterior derivative of (|1.43|) . 

d@ a b - © Q / A r f b + T a f A / 6 = 0, 

that is 

® a b[cd\e] =0, — > R a b[cd;e] = 0- 

This equation says nothing about the Riemann tensor, R a bcd directly. 
It says that certain combinations of the derivatives of the expressions 
on the right hand side of ( (|1.41|) ) are linear combinations of these same 
expressions. 

©ab[cd|e] + D S [ c d® e]sab — ^ & a[c®de]sb ~ ^ b[c® 'de]as = 0. (1-45) 

These are the true Bianchi identities. A consequence of this is that 
if the components of the Riemann tensor are thought of as variables, 
along with the components of the metric and the base forms, then these 
variables have to satisfy 

Rab[cd\e] = -^Rbe[c T df]- ( L46 ) 
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